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Abstract 



►^ , The photon polarization operator in superstrong magnetic fields 

ly^ ' induces the dynamical photon "mass" which leads to screening of 

0^ ■ Coulomb potential at small distances z <^ 1/m, m is the mass of an 

electron. We demonstrate that this behaviour is qualitatively different 
from the case of D = 2 QED, where the same formula for a polariza- 
Q ■ tion operator leads to screening at large distances as well. Because of 

^^ , screening the ground state energy of the hydrogen atom at the mag- 

netic fields B ^ rr? je^ has the finite value Eq = — me^/21n^(l/e^). 
>. 

i-3 1 Introduction 



The Larinour radius of the electron orbit au = l/veB is much smaller than 
Bohr atomic radius a^ = l/(me^) for homogenius magnetic fields B ^ m^e^ 
(we are using Gauss system of units, where e"^ = a = 1/137; also in all 
formulas h = c = 1). It is natural to look for the atomic energy levels in 
such strong magnetic fields studing the influence of Coulomb potential on the 
electrons occupying Landau levels |1J. A strong magnetic field confines an 
electron in the transverse direction while in the longitudinal direction an elec- 
tron is bound by the weak Coulomb field of a nucleus. The cigar-shape wave 
function of an electron is formed with transverse size which equals Larmour 
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radius and longitudinal size which is by ln(a^/a|^) = ln{B/m'^e^) smaller 
than Bohr radius. The ground state energy is larger than Rydberg constant 
by the square of the same logarithm: Eq = — (me^/2)ln {B/m?e^). One 
can easily get this logarithmic factor from the fact that in one-dimensional 
Coulomb potential energy diverges logarithmically at small distances. The 
divergency is regularized at the longitudinal distances which equals au, where 
one-dimensional motion converts to a three-dimensional one. Atomic levels 
in such strong magnetic fields were found numerically in [2] (see also [31 H]). 

Our purpose is to understand the behaviour of the energy levels with the 
growth of a magnetic field. The point is that at superstrong magnetic fields 
B > m? /e^ the polarization operator insertions into the photon propagator 
induce the dynamic photon "mass" m^ ^ e^B [SI E] • One would expect that 
the photon mass should screen Coulomb potential and shift energies of the 
atomic levels found in tree approximation. 

Dirac equation spectrum in a constant homogenious magnetic field looks 
like [7]: 

el = m'^ + pI + {2n + l)eB + aeB , (1) 

where n = 0,1,2,..., a = ±1 and the field is directed along axis zu In 
the magnetic fields we are interested in £„ > m/e, and electrons are ultra- 
relativistic. The only exception is the lowest Landau level (LLL) which has 
n = 0, a = —1. The energy of LLL electron equals its mass for Pz = and the 
consideration of the nonrelativistic electron motion along z axis is selfconsis- 
tent. LLL is interesting both practically and theoretically. An analog of the 
critical electric field Ecr = m'^/e is the magnetic field Bq = th? /e = 4.4 ■ 10^^ 
gauss. Two orders larger superstrong fields B > m? /e^ can exist at special 
neutron stars named magnetars. The temperature of an outer magnetar layer 
is not enough to populate the excited Landau levels and one can observe the 
transitions among the states to which LLL is splitted at the electric field of 
the nucleus. Freezing of the ground state energy in the superstrong magnetic 
fields discussed in the paper leads to the upper bound on the spectra of pho- 
tons radiated from magnetarso To study the stability of the huge magnetic 
fields [9] one should also know the energy of the ground state as a function 
of a field. 

So we are studying the energies of the states to which LLL splits in the 



^Tliis spectrum with the substitution of 2n + 1 + u by 2j, j — 0, 1, 2, ... was found by 
I.I. Raby [8,. 

■^I am grateful to S.I. Blinnikov for the discussions of magnetar physics. 



presence of an atomic nucleus. 

Since the electron at LLL moves along z axis we will study in section 2 
QED at D = 2: the behaviour of electrons in two-dimensional space-time. 
The coupling constant g has dimension of mass, so Coulomb potential as 
a function oi\z\ depends on two dimensionfuU parameters: g and electron 
mass m. We will obtain the approximate analytical formula for Coulomb 
potential m. d = D — \ = \ which takes into account the photon polarization 
operator. We will see that for large g (or small vti) g ^ m Coulomb potential 
is screened. 

In section 3 we will consider the physical case, D = 4 QED. The analog 
of the coupling constant squared g^ in the real world is the product e^B. The 
polarization operator in the magnetic fields B ^ m'^/e at fcif(= k'^) ^ eB 
practically coincides with the one obtained in section 2 [12]. Nevertheless 
the screening at large distances \z\ S> 1/m does not occur: at \z\ 3> l/m 
we get a purely Coulomb potential $(z) = e/\z\. The screening occurs at 
small distances, and its influence on the ground state energy is determined 
in section 4. The results similar to those presented in sections 3 and 4 were 
obtained in [TU] with the help of the numerical calculations. 

2 D = 2 QED: screening 

Coulomb potential in the coordinate space could be obtained by the Fourier 
transformation of the 00-component of the photon propagator in momentum 
representation at momentum k^ = (0, k\\). We designate the space-like com- 
ponent of momentum by /cy, which will be natural in the case d = 3, see 
below. The series of Feynman diagrams for the photon proparator is shown 
in Fig. 1 which corresponds to the following equations: 

$(fc) = Aoik) = ^; ^ = Ao = Doo + DooUooDoo + ... 
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Fig 1. Modification of the Coulomb potential due to the dressing of the 

photon propagator. 



Summing the series we get: 
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P + n(A;2 



T-r _ / _ kflf^U 



U{k' 



(2) 



where U^i, is the photon polarization operator at one loop. Instead of cal- 
culating the fermion loop we can take an expression for 11 obtained in the 
dimensional regularization method |TT], substitute Z) = 2 in it and divide it 
by two, because in two dimensions the traces of 7-matrices are proportional 
to 2 instead of 4: 



n(P) = V 
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t(i + t) 



IniVT+t + Vi) - 1 



= -^g'P{t) 



(3) 



t = -k^/Am^ - 

- a well-known result; for example see [12]. |j 

Taking fc = (0, A;||), A;^ = — /cj? for the Coulomb potential in the coordinate 
representation we get: 



Mz) 



Ang 



^ikuz 



dkii/2n 



k?,+Ag-^P{k?JAm^) ' 

and the potential energy for the charges +g and —g is finally: 

V{z) = ~g^{z) . 



(4) 



(5) 



The calculation of ^{z) would be simplified if we were interested in the 
correction to the potential ~ g^. Expanding denominator of (^ and taking 
into account the first two terms we would deform the integration contour 
in the plane of complex ku in such a way, that the integration result will 



•^It was demonstrated in [12| that in strong magnetic fields photon polarization operator 
is dominated by the electron from LLL and is essentially given by the D — 2 expression. 



be given by the residue at k\\ = and integration of discontinuity of -P(fc|) 
which equals the imaginary part of it. This is analogous to what is done in 
the textbook [13] when the Uehling-Serber correction to Coulomb potential 
in c? = 3 is calculated. However to obtain the photon mass we should take 
into account the whole infinite series - mass is not generated in a finite order 
of the perturbation theory. Discontinuity of the integrand of (jl]) is not equal 
to ImP and the simplification of the integration does not occur. 
Asymptotics of P(t) are: 



P{t) 






t< 1 
t> 1 



(6) 



Let us take as an interpolating formula for P{t) the following expression: 

2t 



Pit) 



(7) 



3 + 2t 

We have checked that the accuracy of this approximation is not worse 
than 10% for the whole interval of t variation, < t < oo. Substituting(I7j) 
in (jlj) we easily perform the integration: 
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2g^l'im^ 



kl + 6m2 + 4(yf2 
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(8) 



^V3m2 



=exp(— A/Gm^ + 4(7^1^1) 



In the case of heavy fermions ivn ^ g) the potential is given by the tree 
level expression; the corrections are suppressed as g^ jvrP': 

*W|„», =-2-9n(i + o(^)) . (9) 

In case of light fermions [m <^ g) the second term which describes Yukawa 
potential m.d = \ dominates at the distances \z\ < {l/g)\n{g/m). At larger 
distances the first term dominates; a coupling constant is suppressed by the 
factor 3m'^/2g'^ with respect to the tree level expression: 
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(10) 



The dependence of the potential energy of the two opposite charges (^ 
on the distance between them is shown in Figure 2o 



Fig. 2. The potential energy of two opposite charges in D = 2 QED in the 
case g ^ ni. The dashed line shows V{z) for g <^m. 

In case m = a hnear term disappears and the potential is determined by 
the photon with the mass m-y = 2g exchange (the Schwinger model: D = 2 
QED with massless fermions [H]). For massive fermions at the distances 
larger than ln[g/m)/g we obtain a linear potential with the coupling constant 
diminished by the factor 3m'^/2g'^ . 

3 D = 4: QED in a strong magnetic field: screen- 
ing at 1 2; I < 1/m 

In order to derive the potential of the pointlike charge in the realistic case of 
D = A QED with d = 3 space-like dimensions in the external magnetic field 
we should know the expression for a polarization operator in this field. There 
are many papers where the polarization operator in the external homogeneous 
field was calculated, see [15] - [17]. The expression for the polarization 
operator radically simplifies in the magnetic field which is so strong, that the 
Landau levels spacing is considerably larger than the electron mass, B ^ 
Bq = ni? /e and at the longitudinal (parallel to the magnetic field) photon 
momentum fc| < eB, see Eqs.(1.19), (1.22) and (5.2) in [IT]. With the help 
of these formulas we get: 



k^ + X2{k^) 



^I am grateful to A.V. Smilga who noted privately that in the case of light fermions in 
D = 2 QED a massive pole in a photon propagator emerges. 
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Aire 



(A:f + A;i)(l + tln(^)) + -i^exp 



_J±\p(J± 

2eB / I 4m2 



(11) 



A;^, + A;^, kz = k\\ and the magnetic field is 



where k = {0,kx,ky,kz), kj_ 
directed parallel to the z axis. This formula is very similar to the analogous 
formulas from the previous section, see Eqs. ([2]), (|3]). The difference is in an 
extra small term ~ aln(ei?/r?7,^), which we will not take into account in the 

fc2 



future consideration, and in the factor 2e B/nexp i —3^ ) which substitutes 

Ag'^. The dependence of function P on kj is the same as that in the case of 
D = 2 QED. There is also an extra term k^ in the denominator and to obtain 
the potential in the coordinate representation we should integrate over k± as 
well. 

With the help of interpolating formula P{t) from section 2 we obtain: 

e'''\\'dk\\d^k±/{2nf 



Hz) 



Ane 



k^ + k^: + 



2e-^B 



exp( 



kl/{2eB)){kl/2m?)/{'i + kl/2m?) 



(12) 
where the integration is performed in the cylindrical coordinates and we are 
looking for the potential along axis z, since it is the potential which bounds 
an electron in the atom. 

We manage to find the asymptotic behaviour of *^{z) for z much larger 
and much smaller than Compton wave length of the electron. For large 
distances \z\ ::^ — in the integral flT^ the region |A;||| ^ m is important and 
for the magnetic field B ^ Bq we get fcj? <^ eB and the expression for P we 
are using is correct. For small \k\\\ we get: 
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\z\ > 
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(13) 



It is convenient to integrate over k\\ closing the integration contour in an 

upper (lower) semiplane of the complex k\\ and taking k\\ = ik±/ 
In this way we obtain: 
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(14) 



The integral over k± converges at k± < Ve^B, that is why the residue 
was situated at k^\ ^ m, where the approximate formula for P we used is 
valid. At the mentioned values of k_j_ the exponent inside the square root can 
be substituted by one and finally we obtain: 



^{z) 



z > ^ 
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\z\ 



V{z) 



m 



e 
\z\ 



(15) 
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- the usual Coulomb potential. Strong screening which we obtain in d 
at the distances 1^1 ^ 1/m in a realistic case d = 3 does not occur. 

To find a potential at short distances l^;] -C ^ let us substitute m = in 
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eB 



-k±dk± .(16) 



Calculating the potential at \z\ ^ XjyeB we observe that the integral over 
fc_L is determined by the integrand at k^ <^ a/^B. So we took residue at 
k\\ ^ k± <^ \/eB and the approximate expression for P was used in the 
domain where it is valid. Performing integration over k± we get: 



^(z) 



:^ >-2 > ^ 



eB 



V{z) 




k\dk\ 



(17) 



At the distances which are smaller than Compton wave length we obtain 
screening of the potential which corresponds to the photon mass m? = 
2e^B/7i. Coulomb potential is screened for the superstrong magnetic fields 
B > m^/e^. 



4 The energy of the ground state of the hy- 
drogen atom in the superstrong magnetic 
fields B > ml/e' 
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According to papers [21 E] in the magnetic fields B > rri^e^ the ground state 
energy of the hydrogen atoms equals Eq = — (me^/2)ln {B/m?e^) and at 
B = rn^/e^ it equals Ecr = — (me^/2) In (1/e^). For larger magnetic fields 
the screening of the Coulomb potential at the distances |z| < — occurs. Let 
us demonstrate that the screening leads to the freezing of the energy - it 
does not diminish with the growth of the magnetic field. 

To find the ground state energy we use the following equation [3]: 



En = -2m 



We split the integral into two parts: from 1/m to as, where the screening is 
absent (large z), 

as 2 

h = - I -dz = -e^ In (l/e^) (19) 

l/m 




and from the Larmour radius an = l/^feB to l/m, where the screening 
occurs (small z): 

l/m 2 

/2 = - j — exp(-V^z)rfz = -e^ln(l/e) . (20) 

l/v^ 

Finally we get: 

Eq = -meV21n2(l/e^) = -me^/2 x 220 (21) 

and the energy level freezing occurs. The numerical estimates of Shabad and 
Usov give 73.8 x 4 f» 295 instead of 220, see Eq. (14) in Phys. Rev. Lett, 
paper pH] . 

When B increases further Larmour radius approaches the size of a proton. 
This happens at 1/y/eB Ri l/rrip, m^ = 770 MeV, B ^ 10^° gauss. Taking 
into account the proton formfactor we get that for larger fields I2 does not 



contribute to the energy, factor 220 in (!?!]) should be substituted by 100: the 
ground level goes up. 

Without screening / = — e^ \io.{aB / clh) , Eq = — (me^/2) In (B/m^e^) as it 
was stated in the beginning of this section. 

5 Conclusions 

The photon polarization operator leads to modifications of the atomic en- 
ergy levels. The famous example is its contribution to the Lamb shift, the 
difference of the energies of 2si/2 and 2pi/2 levels of hydrogen. They are 
numerically small loop corrections to the values of energies determined by 
the tree level potential. The role of the photon polarization diagram in the 
superstrong magnetic fields B > w? /e^ = 6.2 ■ 10^^ gauss is qualitatively dif- 
ferent. It determines the behaviour of the ground state energy: the formula 
obtained at tree level becomes invalid and the growth of the coupling energy 
with B terminates at -B ~ m? /e^. Screening of Coulomb potential should be 
more important for the energies of even excited states which are more sensi- 
tive to the shape of the potential at small distances [18] . Degeneracy of even 
and odd excited states in the limit B =^ oo is not lifted by the screening. 
We study the analogy of the electric potential in rf = 1 QED with massive 
electrons and in c? = 3 QED in strong magnetic fields B > Bq = m? /e 
which originates from the coincidence of the polarization operators in these 
cases. A simple analytical expression which equals the polarization operator 
with 10% accuracy enables us to obtain an approximate formula for the elec- 
tric potential of the point charge m. d = 1 QED with massive fermions and 
asymptotics of the potential in c? = 3 QED. In d = 1 QED for a coupling 
constant g larger than a fermion mass m a tree level formula is modified at 
\z\ > 1/g. In d = 3 QED a tree level formula is modified at the distances 
1/rra > \z\ > 1/Ve^B while at large distances \z\ > 1/m Coulomb law is 
valid. 

Analogous results for D = A were obtained in |10j . 

The other aspect of the Coulomb potential in the strong magnetic field is 
investigated in paper [19]: it is supposed that fermions obtained their mass 
due to a magnetic field (dynamical fermion mass). 

I am grateful to S.I. Blinnikov, V.A. Novikov, L.B. Okun, V.S. Popov, 
and A.V. Smilga for useful discussions and to A.I.Rez, who brought to my 
attention Phys. Rev. publication [TO] . 
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